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HYBRID SPACE-TIME DIVERSITY BEAM FORMING SYSTEM 



Field of the Invention 

The present invention relates to hybrid space-time diversity beam forming, 
and is particularly concerned with applications to mobile wireless systems. 
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Background of the Invention 

A GSM network is composed of several functional entities, whose functions 
and interfaces are defined. The GSM network can be divided into three broad parts. 
The Mobile Station is carried by the subscriber; the Base Station Subsystem controls 
the radio link with the Mobile Station. The Network Subsystem, the main part of 
which is the Mobile services Switching Center, performs the switching of calls 
between the mobile and other fixed or mobile network users, as well as management 
of mobile services, such as authentication. The Mobile Station and the Base Station 
Subsystem communicate across the Urn interface, also known .as the air interface or 
radio link. The Base Station Subsystem communicates with the Mobile service 
Switching Center across the A interface. The following table briefly describes the 
GSM development history 
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Table 1 GSM Development History 





Phase 1 


Phase 2 


L^hase 2+ 


Specifications 
development 


1987-1991 


1992-1995 


1995-1999 


Trial period 


1992-1998 


1998-1999 


1998-2000 


Voice transmission 
modes 


Full Rate 


Full rate, half rate 


Full rate, half rate, 
enhanced full rate 


Data transmission 


Full rate (9.6kbps max) 
Half rate (4.8kbps max) 


Full rate (9.6kbps 
max) 

Half rate (4.9kbps 
max) 


Full rate (1 1.4kbps max), 
Half rate (4.8kbps max), 
Packet transmission 
(T6.8kbps max), 
High-speed circuit 
exchange (T8.8kbps 
max) 
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Supplementary services 


Call transfer, call 
limiting 


Call transfer, call 
limiting, hold, caller 
ID, service of 
charge, conference 
calls (6 parties max), 
etc. 


Enhanced versions of 
Phase 2 services 


Short message service 


Bidirectional text 
transmission, cell- 
specific announcements 


Category-specific 
text transmission 
(basic, voice mail 
notification 


Phase 2 capabilities 
plus data packet 
transfer and two-byte 
code transfer 


Other 


Global roaming 5V BIM 
card 


Global roaming, 3- 
SV SIM card, 
DCS 1800 
specifications 


Global and domestic 
roaming, SIM tool kit, 
multiband GSM, DCS 
support, 

interoperability with 
DECT 



GSM Channel structure 

5 GSM is an FDD-TDMA system, each carrier occupies 200 kHz which is time- 

shared by 8 time slots or users. The structure of the most common timeslot burst is 
shown in Fig. 1. The structure 10 includes three tail bits 11, 58 information bits 12, 
26 training bits 13, a second group of 58 information bits 14, three tail bits 15 and 
8.25 guard bits 16. 

10 

A total of 156.25 bits is transmitted in 0.577 milliseconds, giving a gross bit 
rate of 270.833 kbps. There are three other types of burst structure for frame and 
carrier synchronization and frequency correction. The 26 bit training sequence 13 is 
used for slot timing and equalization,, as described below. The 8.25 bit guard time 16 
15 allows for some propagation time delay in the arrival of bursts. 

Referring to Fig. 2, there is illustrated the GSM Frame Structure. Each group 
of eight time slots is called a TDMA frame 20, which is transmitted every 4.615 
ms. TDMA frames are further grouped into multiframes to carry control signals to the 
20 designated slots/blocks. There are two types of multiframe 21 and 22, containing 26 

or 51 TDMA frames respectively. The 26 frame multiframe contains 24 Traffic 
Channels (TCH) and two Slow Associated Control Channels (SACCH) that supervise 
each call in progress. The SACCH in frame 12 contains eight channels, one for each 
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of. the eight connections carried by the TCHs. The SACCH in frame 25 is not 
currently used, but will carry eight additional SACCH channels when half rate traffic 
is implemented. A Fast Associated Control Channel (FACCH) works by stealing 
slots from a traffic channel to transmit power control and handoversignalling 
5 messages. The channel stealing is done by setting one of the control bits in the time 

slot burst to indicate the slot carries traffic data or signaling message. 

In addition to the Associated Control Channels, there are several other control 
channels, which except -for the Standalone Dedicated Control. Channel are 
10 implemented with broadcasting channel (say, in time slot 0) of specified TDMA 

frames in a 51 frame multiframe 22, implemented on a no hopping carrier frequency 
in each cell. The control channels include: 

• Broadcast Control Channel (BCCH): Continually broadcasts, on the downlink, 
information including base station identity, frequency allocations, and frequency 
hopping sequences etc. 

• Standalone Dedicated Control. Channel (SDCCH): Used for registration, 
authentication, call setup, and location updating. Implemented on a time slot, 
together with its SACCH, selected by the system operator. 

• Common Control Channel (CCCH): Comprised of three control channels used 
during call origination and call paging. 

• Random Access Channel (RACH): A slotted Aloha channel to request 
access to the network 

• Paging Channel (PCH): Used to alert the mobile station of incoming call. 

• Access Grant Channel (AGCH): Used to allocate an SDCCH to a mobile 
for signaling, following a request on the RACH. 

GSM Frame Structure and Numbering 

The TDMA technique means that the data are multiplexed in time blocks, in 
15 GSM one uses slots or frames. The frames are then grouped into multi-traffic frames 

(26 normal frames) or multi-control frames (51 normal frames). These two types of 
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multi-frames are grouped again into super-frames 23 and 24 with 26x51 normal 
frames. So the numbering of each frame is start from 0 to 2048x26x51 - 1 = 
2,715,647. The hyperframe 25 is the largest cycle and is repeated in the network. 

Mobile needs to inform the network the signal strength information from 
surrounding cells by signaling channels such as Associated Control Channel. 

Bandwidth required to send the measurement information is 1/24 of that 
required for voice. 



The following table summarizes the structure and numbering. 





Traffic Channel 
(TCH) 


Consists of: 


Control 

Channel (CCH) 


Consists of: 


TDMA 
frame 


4,6 ms ' 


8 Timeslots 


4.6 ms 


8 Timeslots 


Multiframe 


120 ms 


26 TDMA 
frames 


234 ms 


51 TDMA 
frames 


Superframe 


6 s 120 ms 


51 Multiframes 


6 s 120 ms 


26 Multiframes 


Hyperframe 


3 h 28 m 53 s 760 
ms 


2048 

Superframes 


3 h 28 m 53 s 
760 ms 


2048 

Superframes 



Frame Synchronized Transmission 

For each carrier with 200kHz spectrum, a maximum of 8 users may share, this 
spectrum in time domain. The biggest issue in uplink is to guarantee that those signals 
from randomly appearing users do not overlap each other in time slots while keeping 
the overhead small. GSM has designed a so called 'Timing Advance' to resolve this 
issue. Referring to Fig. 3 there is graphically illustrated the timing advance zones of 
GSM. 
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In GSM, the targeting range is up to 35 km radius 26 of the Base Station (BS). 
A round trip delay for the far most point is 233.3 jis. In order to guarantee uplink 
signals of the same carrier will not overlap each other, a guard time is necessary for 
each mobile. To handle 233.3 jis round trip delay, we need 252 ins (or 68.25 bits) 
5 which is bad in terms of the spectrum efficiency. So only the RACH has this luxury. 

The traffic channels cannot waste that much bandwidth. The 'timing advance' used in 
GSM causes a mobile to transmit earlier than just 3 slot delay relative to downlink 
timing. 

10 Referring to Fig. 4, there is graphically illustrated the GSM timing advance 

processing. The BS transmits a burst 30. The MS receives it 32 T(l) later. The MS 
transmits a burst 34, 3 slots later. The BS receive it 36 T(2) later. When a base 
station detects a RACH, it meanwhile measures the round trip delay T(l) +T(2) and 
derives the delay from mobile to BS. This delay is called timing advance and is 

15 signaled to mobile to adjust its transmitting time 38 to 3 slot minus that timing 

advance. In summary, the maximum round trip delay 233.3 \ls is quantized to a 6 bit 
number, so 64 steps (0-63) possible timing advance! Each step advances the Timing 
by one bit duration 48/13 (233.3/63) jis which is about 70/64 = 1,0938 km. This 
timing advance process can reduce the guard time to be 8.25 x 3.7 = 30.525 |is rather 

20 than 252 |is required for the maximum round trip. 

64 steps allows for compensation over a maximum one way propagation time 
of 3 1.5 bit periods i.e. 116.3 \is (i.e. a maximum distance of - 35 km) 

25 • Initial Timing advance: BS instructs the MS who advances its burst 
transmission by a time corresponding to round trip delay. The maximum timing 
advance value is 63. (GSM 03.30 defines how PLMN deals with MS when the timing 
advance value is greater then 63). 

30 * Tracking Mode Timing Advance: The BS continuously monitors the delay of 

the normal bursts sent by MS. If the delay changes by more than 1 bit period, the 
timing advance shall be advanced or retarded only 1 bit and the new value signaled to 
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MS. The purpose of restricting the timing advance to 1-bit period each time is to 
simplify the implementation in BS. However, BS may use "large" stepsize (ref GSM 
05.10> 

5 The Timing Advance is used to compensate for the time it takes a RF signal to 

go at the speed of light between the BTS and MS. The maximum BTS radius of 35 
km is divided into 64 TA steps (This means 547 meters / TA step - As a simplification 
550 meters is used). The TA multiplied with 550 meters gives the minimum distance 
to the BTS. The maximum distance is 550mx (TA +1). A TA value places a BTS in a 

10 circular band 550 meters wide, with an inner radius of (TA x 550) meters. 



Timing 
Advance 


0 


1 


2 


3 


4 


5 




63 


Distance to 


< 


550m- 


1100m- 


1650m- 


2200m- 


2750m- 




35 




550 
















BTS 




1100 m 


1650 m 


2200 m 


2750 m 


3300 m 




Km 




m 

















Speech coding and Channel Coding 

Referring to Fig. 5, there is illustrated the speech coding and channel coding 
1 5 for GSM. GSM is a digital system, so speech signals, which axe inherently analog, 

have to be digitized. The method employed by ISDN, and by current telephone 
systems for multiplexing voice lines over high-speed trunks and optical fiber lines, is 
Pulse Coded Modulation (PCM). The output stream from PCM is 64 kbps, too high a 
rate to be feasible over a radio link. The 64 kbps signal contains much redundancy, 
20 although it is simple to implement. The GSM group studied several voice coding 

algorithms on the basis of subjective speech quality and complexity (which is related 
to cost, processing delay, and power consumption once implemented) before arriving 
at the choice of a Regular Pulse Excited - Linear Predictive Coder (RPELPC) with a 
Long Term Predictor loop. Basically, information from previous samples, which does 
25 not change very quickly, is used to predict the current sample. The coefficients of the 

linear combination of the previous samples, plus an encoded form of the residual, the 
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difference between the predicted and actual sample, represent the signal. Speech is 
divided into 20 millisecond samples, each of which is encoded as 260 bits 40, giving a 
total bit rate of 13 kbps. 

5 Due to natural or manmade electromagnetic interference, the encoded speech 

or data transmitted over the radio interface must be protected as much as is practical. 
The GSM system uses convolutional encoding and block interleaving to achieve this 
protection. The exact algorithms used differ for speech and for different data rates. 
The method used for speech blocks is described below. 

10 

Recall that the speech codec produces a 260-bit block for every 20 ms speech 
sample. From subjective testing, it was found that some bits of this block were more 
important for perceived speech quality than others. The bits are thus divided into 
three classes: 

15 

Class la 50 bits - most sensitive to bit errors 
Class lb 132 bits - moderately sensitive to bit errors 
Class II 78 bits - least sensitive to bit errors 

20 Class la bits have a 3 bit Cyclic Redundancy Code added for error detection. 

If an error is detected, the frame is judged too damaged to be comprehensible and it is 
discarded. It is replaced by a slightly attenuated version of the previous correctly 
received frame. These 53 bits, together with the 132 Class lb bits and a 4-bit tail 
sequence (a total of 189 bits), are input into a 1/2 rate convolutional encoder of 

25 constraint length 4. Each input bit is encoded as two output bits, based on a 

combination of the previous 4 input bits. The convolutional encoder- thus outputs 378 
bits, to which are added the 78 remaining Class EL bits, which are unprotected. Thus 
every 20 ms speech sample is encoded as 456 bits 42, giving a bit rate of 22.8 kbps. 
To further protect against the burst errors common to the radio interface, each sample 

30 is diagonally interleaved. The 456 bits output by the convolutional encoder are 

divided into 8 blocks of 57 bits 44, and these blocks are transmitted in eight 
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consecutive timeslot bursts 46. Since each timeslot burst can cany two 57-bit blocks, 
each burst carries traffic from two different speech samples. 

Recall that each timeslot burst is transmitted at a gross bit rate of 270.833 
kbps. This digital signal is modulated onto the analog carrier frequency, which has a 
bandwidth of 200 kHz, using Gaussian filtered Minimum Shift Keying (GMSK). 
GMSK was selected over other modulation schemes as a compromise between 
spectral efficiency, complexity of the transmitter, and limited spurious emissions. 
The complexity of the transmitter is related to power consumption, which should be 
minimized for the mobile station. The spurious radio emissions, outside of the 
allotted bandwidth, must be strictly controlled so as to limit adjacent channel 
interference, and allow for the coexistence of GSM and the older analog systems (at 
least for the time being). 

Referring to Fig. 6, there is illustrated in a block diagram a GSM transmitter. 
The GSM transmitter main 50 includes a voice Codec 52 coupled to a voice digityzer 
54, and outputting to a channel encoder and interleaver 56. The output of which is 
input to MUX 58 then GMSK modulator 60 and channelizer 62 for transmission. 

20 Modulating symbol rate 

The modulating symbol rate is 1/T = 1 625/6 ksymb/s (i.e. approximately 
270.833 ksymb/s), which corresponds to 1 625/6 kbit/s (i.e. 270.833 kbit/s), 
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Start and stop of the burst 

Before the first bit of the bursts as defined in GSM 05.02 [3] enters the 
modulator, the modulator has an internal state as if a modulating bit stream consisting 
of consecutive ones (d t = /) had entered the differential encoder. Also after the last bit 
of the time slot, the modulator has an internal state as if a modulating bit stream 
consisting of consecutive ones (d f = J) had continued to enter the differential encoder. 
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These bits are called dummy bits and define the start and the stop of the active and the 
useful part of the burst as illustrated in figure 1. Nothing is specified about the actual 
phase of the modulator output signal outside the useful part of the burst. 

Differential encoding 

Each data value <// =[0,1] is differentially encoded. The output of the 
differential encoder is: 

dfdtQdt-i (^e {0,1}) 

where © denotes modulo 2 addition. 

The modulating data value bi input to the modulator is: 

b i =l-2d l (6 <e {-l,+ l}) 

Training Sequence 

There are eight training sequences are designed in GSM. system. Each slot/user 
has a unique training sequence as a midamble. Each training sequence has 26 bits and 
the sequences were designed to have a good cross correlationship. 

Slot-1 Training = [0,0,1,0,0,1,0,1,1,1,0,0,0,0,1,0,0,0,1,0,0,1,0,1,1,1]; 

Slot-2 Training = [0,0,1,0,1,1,0,1,1,1,0,1,1,1,1,0,0,0,1,0,1,1,0,1,1,1]; 

Slot-3 Training = [0,1,0,0,0,0,1,1,1,0,1,1,1,0,1,0,0,1,0,0,0,0,1,1,1,0]; 

Slot-4 Training = [0,1,0,0,0,1,1,1,1,0,1,1,0,1,0,0,0,1,0,0,0,1,1,1,1,0]; 

Slot-5 Training - [0,0,0,1,1,0,1,0,1,1,1,0,0,1,0,0,0,0,0,1,1,0,1,0,1,1]; 

Slot-6 Training = [0,1,0,0,1,1,1,0,1,0,1,1,0,0,0,0,0,1,0,0,1,1,1,0,1,0]; 

Slot-7 Training = [1,0,1,0,0,1,1,1,1,1,0,1,1,0,0,0,1,0,1,0,0,1,1,1,1,1]; 

Slot-8 Training = [1,1,1,0,1,1,1,1,0,0,0,1,0,0,1,0,1,1,1,0,1,1,1,1,0,0]; 
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Filtering 

The modulating data values as represented by Dirac pulses that excite a 
linear filter with impulse response defined by: 
5 g{t) = Kt)*reet(^ 

where the function red (x) is defined by: 



and * means convolution. h(t) is defined by 

exp 

10 A(0 = 



recti ~- = 0 otherwise 



f ^ 



where a«tfJ?r«03 

where B is the 3 dB bandwidth of the filter with impulse response h(t) y and T is the 
duration of one input data bit. This theoretical filter is associated with tolerances 
1 5 defined in GSM 05.05 [4]. 

Output phase 

The phase of the modulated signal is: 

f'-iT 
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where the modulating index h is 1/2 (maximum phase change in radians is n/2 per 
data interval). 

The time reference V = 0 is the start of the active part of the burst as shown in figure 
1 . This is also the start of the bit period of bit number 0 (the first tail bit) as defined in 
GSM 05.02 [2]; 



Modulation 

The modulated RF carrier, except for start and stop of the TDMA burst may 
therefore be expressed as: 

10 *T/n-J*S 



where E c is the energy per modulating bit, f 0 is the center frequency and <p 0 is a 
random phase and is constant during one burst. 

GSM Receiver Chain 

15 Referring to Fig. 8, there is illustrated in a block diagram a GSM receiver. 

The GSM receiver 70 includes a splitter 72, a demodulation 74, a DEMUX 76 a 
channel decoder and deinterleaver 78 the output of which is coupled to a speech 
decoder and digital to analogue converter (DAC) 80. At the 900 MHz range, radio 
waves bounce off everything - buildings, hills, cars, airplanes, etc. Thus many 

20 reflected signals, each with a different phase, could reach an antenna. Equalization is 

used to extract the desired signal from the unwanted reflections. Equalization works 
by finding out how a known transmitted signal is modified by multipath fading, and 
constructing an inverse filter to extract the rest of the desired signal. This known 
signal is the 26bit training sequence transmitted in the middle of every time slot 

25 burst. The actual implementation of the equalizer is not specified in the GSM 

specifications. 
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Matching Filter 

After down converter and ADC, the received sampled data is usually flittered 
by an anti-alising filter, which will decide the sample phase so that to decimate to 
symbol rate data. > 

5 

Time Synchronization 

The 26 known bits are used to do a correlation base search and to find the slot 
boundary 



10 Channel Estimation 

Again, using the known 26 bits to do the channel estimation. A) Correlation 
based, B) Least mean square error solution. 



Viterbi Equalization 

With the estimated channel (usually 5 Or 7 taps), enumeration search on the 
possibility set is pursued to find the possible bit sequence, which will be the input for 
channel decoder. 

GSM has become a very successful wireless scheme with repaid subscriber 
growth in its areas of deployment. In fact growth has been in such a rapid rate that 
the capital costs of installed infrastructure has not had enough time to pay for itself 
before further equipment deployment must be made in order to service now 
subscribers. This growth related problem has been most acutely felt in urban areas. 
Consequently, wireless service providers need to find a way to provide services to a 
growing customer base without growing the infrastructure at the same rate. 
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Summary of the Invention 

An object of the present invention is to provide a hybrid space-time diversity 
beam forming system. 

la accordance with an aspect of the present invention there is provided a 
method of beam forming comprising the steps of: in an appliquS intelligent antenna 
system, monitoring broadcast channels of a mobile wireless base station; monitoring a 
frequency burst broadcast by the base station and synchronizing the applique system 
in frequency; monitoring a synchronization burst in the broadcasting channel and 
synchronizing the applique system with the mobile wireless base station in time. 

Conveniently, a step of determining an angle of arrival includes the step of 
determining a co variance matrix XX, where X is given by: 

£1(62) £2(62) £3(62) £4(62)" 
£1(63) £2(63) £3(63) £4(63) 
£1(64) £2(64) £3(64) £4(64) 

£1(87) £2(87) £3(87) £4(87) 



AndXX = X*-X 

1 5 Conveniently, the step of determining the angle of arrival includes the step 

of forming a Hermitian Toeplitz matrix by using XX with the following procedures 

Zo - [XX(1,1) + XpC(2,2) + XX(3,3) + XX(4,4)] / 4; 
Zi = [XX(1 ,2) + XX(2,3) + XX(3,4)] / 3; . 
Z2 = [XX(l,3)+XX(2,4)]/2; 
20 Z 3 = XX(1,4) 



An advantage of the present invention is increasing the subscriber capacity of 
an existing base station. 
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Brief Description of the Drawings 

The present invention will be further understood from the following detailed 
description with reference to the drawings in which: 

Fig, 1 illustrates a known CSM common burst slot structure; 

Fig. 2 illustrates a known GSM frame structure; 

Fig. 3 graphically illustrates known timing advance zones for GSM; 

Fig. 4 graphically illustrates known timing advance processing for GSM; 

Fig. 5 illustrates known speech coding and channel coding for GSM; 

Fig. 6 illustrates in a block diagram a known GSM transmitter chain; 

Fig. 7 graphically illustrates a known relationship between active part of burst, 
tail bits and dummy bits; 

Fig. 8 illustrates in a block diagram a known GSM receiver chain; 

Fig. 9 illustrates a four-element linear antenna array system in accordance 
with an embodiment of the present invention; 

Figs. 10 and 11 illustrate a simplified receiver portion and transmitter portion, 
respectively of an appliqud intelligent antenna system for use with the antenna array 
system of Fig. 9 in accordance with an embodiment of the present invention; 

Fig. 12 illustrates a receiver for the antenna array system of Fig.. 9 in 
accordance with an embodiment of the present invention; 
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Fig. 13 illustrates in a block diagram a transmitter for the intelligent antenna 
system in accordance with an embodiment of the present invention; 

Figs. 14 and 15 illustrate the detailed arrangement of information in the SCH 
message; and 

5 Fig* 16 graphically illustrates downlink beam patterns; 

Detailed Description of the Preferred Embodiment 

Referring to Fig. 9, there is illustrated a four-element linear antenna array 
10 system in accordance with an embodiment of the present invention. The four element 

linear antenna array system 100 includes antennas 102 spaced a distance d=0.618 
wavelength apart. This is actually the golden number, which we found to be the best 
compromise between traditional linear array antennas and diversity antennas. With 
this golden separation, , we can collect antennae diversity while keeping the merits of 
15 the linear array processing gains. In an embodiment of the present invention the 

antenna element spacing is (5 Yl -l)/2 times the wavelength. 

The four-element linear array system of Fig. 9 is designed to improve the link 
quality of both uplink and down link. The four-element linear array system including 

20 beam-forming technology is implemented for both uplink and downlink, which can be 

integrated into the customer's base station TRU. The principle of a linear array can 
be illustrated in the following array model. Suppose a signal, s(t) = B(t) exp (j27tft) 
with a base band signal B(t), being transmitted from a mobile phone, is impinging the 
linear array. Regarding to the physical array, the signal wave can be regarded as a 

25 plane wave. 

Suppose the distance between the two adjacent elements is d meters, as the 
signal propagate with light speed c, so the signal arrives the 2 nd element will be 
delayed by d x sin (0)/c seconds. Similarly, the signal arrives at 3 rd array element will 
30 be delayed by 2 x d x sin (Q)/c and arrives at 4 th element will be delayed by 3 x d x 



WO 2004/047333 



PCT/CA2003/001747 



16 



10 



sin (6)/c. As the symbol duration (48/13 micro seconds) is much larger then the 
duration of the signal sweeps of the array, the narrow band signal received by each 
element can be regarded as unchanged except its phase. Then we may model the array 
output is 



>,(0" 




1 


x 2 (t) 




expO'2^* x d x sin(0) / c) 


x 3 (f) 


exp(J2tf x 2 x d x sin(0) / c) 


_x 4 (t) 




exp(J2nf x 3 x d x sin(0) / c) 



= *(0 



1 • 

expO'2^r xdx sin(0) / ;i) 
expO'2^ - x 2 x sin(0) / /I) 
exp(J2n x 3 cf x sin(0) / ,1) 

(1) 

where X = c/fis the wavelength. A particular design is to make the distance between 
every two adjacent elements equals to haif wavelength, and then the above equation 
can be simplified as: 



x 2 (t) 



=K0 



1 

exp(y>zrsin(0)) 
exp(y a: x 2 x sin(#)) 
_exp(yVr x 3 x sin(0)) 



(2) 



Note that we have assumed the signal received by the 1 st array element has no 
15 propagation delay. This is only true after the received signal has been synchronized 

with transmission timing. Also we did not count any imperfections such as multipath, 
interference and noise and so on. More generically, we can assume many signals 
along with their multipath copies are simultaneously impinging the array. 

20 For GSM, as it is FDD-TDMA system, only one desired signal and its 

multipath need to be extracted, other impinging signals are interference or noise that 
need to be suppressed as much as possible. We now suppose the desired signal is 
propagated via L paths. After demodulation, the received baseband array output can 
be modeled as 
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1 



exp(y" 2?r x J x sin(#(/)) / X) 
exp(j2n x 2d x sin(6>(/)) / A) 
expC/27ir x 3d x sin(<9(/)) / X). 



+ I(t) 



{*(/)} 
{*(/» 



+ /(*) 



(3) 



where {b (i)} is the transmitted symbol sequence, I(t) is the overall noise and 
interference effect, chk(t) is the multipath channel modulated by the array steering 
vectors. More precisely, 



cA,(0 = £«(/)/»(< -r(0) 
ch z (0 = Z a O ex PO' 2 ^ x d x sin(0(/)) / - 
CA3 (0 = Yi a (0 «p0'2tf x 2d x sin(0(/)) / X)h(t - r(i)) 



ch 4 (0 = 2 «(0 expO*2^ x 3d x sin(0(/)) / - r{l)) 



(4) 



10 



where h (t) is the overall channel impulse response. 



15 



The four-element linear array system fully utilizes the known sequence to 
estimate all the parameters such as slot boundary, angle of arrival (AO A) and time of 
arrival (TO A) of the strongest path, channel impulse response which is used for both 
uplink decoding and downlink transmission. 



20 



Referring to Figs. 10 and 1 1, there are illustrated a simplified receiver portion 
104 and transmitter portion 106, respectively of an applique intelligent antenna 
system for use with the antenna array system of Fig. 9 in accordance with an 
embodiment of the present invention. On the receive portion 104, the antenna array 
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system 100 is coupled to an RF/EF down-ponverter 108 that provides input to a digital 
base-band processing block 110 whose output is applied to an IF/RF up-converter 
1 12. The output of the IF/RF up-converter 1 12 is applied to an existing base station 
114. On the transmit portion 106, the existing base station 1 14 is coupled to an RF/IF 
5 down-converter 116 that provides input to the digital base-band processing block 110 

whose output is applied to an EF/RF up-converter 118. The output of the IF/RF up- 
converter 1 18 is applied to an antenna array system 100. Consequently, both receive 
and transmit portions of the appliqu6 intelligent antenna system are transparent to the 

existing base station. 

» 

Referring to Fig. .12, there is illustrated the receive portion of Fig. 10 in further 
detail. The receive portion 104 includes four analog-to-digital converters (A/D) 120 
each coupled to a respective antenna 102 with outputs coupled to a slot framing and 
timing block 122. The outputs of the slot framing and timing block 122 are each 
applied to a respective mixer 124 and separately applied to a beam weight block 126. 
The outputs of the beam weight block 126 are applied to the mixer 124 whose outputs 
are applied to an adder 128, whose output is applied as input to an IF/RF up-converter 
130. The outputs of framing/timing block 122 are also applied to a four-channel 
estimation block 132 and an angle of arrival (AO A) estimation block 134. The AOA 
estimation block 134 is output to a transmit weight block 136 whose output is coupled 
to a transmitter chain (not shown in Fig. 12). The outputs of the framing and timing 
block 122 corresponding to the outer antennas of the array are coupled to a diversity 
selector 138. In operation, the slot framing/timing block is similar to a GSM timing 
block, i.e. it uses the known 26 bits to do correlation with the received samples and to 
find out the accurate slot boundary. One difference from typical GSM is here we have 
four data flow inputs instead of one, hence the framing/timing here pursues a common 
timing for all four data flows. 

Note that, in this stage, we have made the assumption that the four-element 
antenna array system has acquired the BS timing from the synchronization burst (see 
the Watchdog description herein below for details). 



15 



20 



25 
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The output of this block will be a slot wise data vector with symbol rate or 
over sampling rate. The sampling rate M is specified to be 1, 2 and 4 samples per 
symbol. The complexity for this block is summarized in the following table. 

Table 1 Complexity of Slot Framing 



Searching Window 


Pseudo MPX 


M =2 


M = 4 


7 symbols 


4x7xMx26 = 


1456 


2912 



10 



• Pseudo multiplication means the implementation can be done as an 'adder' 
since the known sequence is a sort of +1, -1, +j, -j. 



The CCIC Beamformer Weights Block 126 is responsible for calculating the 
beamformer weights which are used to combine the four data flows to form the input 
for the known base station receiver 142. 



15 



20 



25 



Space-Time Diversity Beamforming (STDB) Algorithm 

Suppose the framed symbol rate slot data is 



yi 

y^ : \ 

y 3 
y* 



R1(1),R1(2),R1 (3),R1(4) R1(61),RJ(62) > ...,RI(87),R1(88),... ^1(145)^1(146)^1(147)^1(148)^1(149) 

,...,R1(156) 



R2(1)^2)^(3)J^(4),...,R2(61),R2(62),.:,^^ 
...R2(156) 



R3(l)JG(2)JR3(3)Jtt(4),...,R3(61),R3(62)^ 
,...R3(156) 



R4(1),R4(2),R4(3),R4(4) > . . .,R4(61),R4(62),.. .,R4(87)JU(88),. . . JU(145)JR4(14<5),R4(147),R4(148),R4(149) 

...R4(156) 
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Refer to equation (3); the array output at time k can be expressed as 

" cA,(6)Y s(Jk) ^ 

• ch 2 (6) s(k-Y) 

• ch 3 (6) : 







r «*i(Q) 


.^(1) 






ch 2 (0) 


cfc 2 (l) 






ch 3 (0) 








,cA 4 (0) 





where chl, ch2, ch3 and ch4 are the estimated channel impulse response, s(k)'s are 
the transmitted MSK symbols. Especially when we choose those s(k)'s to be the 26 
known training sequence and arrange the array output into a space-time data array as 



y,W y>(k + i) 



y>(k + 25) 

y,(k + 25) 
^(fc + 25)J 



fcft,(0) cA,(l) 

cA,(0) cA,(l) 

cA,(0) C A 3 (1) 

«M0) ch.(l) 



ch,(6)-) 
cA 3 (6) 

<*4(6)A*(*-6) 



s(k-l) 



s(k + 1) 



sQc + 2Sp 
s(k + 24) 



j(A-5) .« *(A + 19), 



+ IN 



10 



Note that in this equation, only IN, the interference plus noise seems unknown 
and we should minimize its affect. So our optimal beamformer will be chosen such 
that . 



15 



min{w T JNJN T w = w T (Y- ChS)(Y - ChS) T w,s.t. || vv || 2 = l}. 

The solution for this minimi2ation problem is again an eigenvalue problem of 
a 4x4 semi-definite positive Hermitian matrix, that has an explicit solution. One way 
to solve this optimization problem is to do an eigen value decomposition for the 4x4 
Hermitian matrix (Y-ChS)(Y-ChS) T . 
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The Complexity of this block is summarized in the following table. 
Table 2. Complexity of Algorithm I 





Complex MPY 


Real MPY 


Real Addition 


Total 

Instructions 


Conj(X) T X 


4*26 = 104 


416 


240 




Inv(XX) 


64 


256 


16 




y 


4*26 = 104 


416 


240 




w 


16 


64 


24 


1672 



Remark: The big difference between the present implementation and others is the 
5 treatment for the correlation matrix XX when it is invertible. 

Beamforming Block 

• Combine the four data vectors by using formula Db(k) = conj(w(l))Rl(k) + 
conj(w(2))R2(k) + conj(w(3))R3(k) + conj(w(4))R4(k), for k = 1, 2, . . 156, 
10 •Up convert to RF at block 140 and theri feed into receiver 142 having down 

connector 144 and receiver block 146. 



The complexity of this block is summarized in the following table. 
Table 3 Complexity of Beamforming Block 





. Complex 
MPY 


Real MPY 


Real Addition 


Total 

Instruction 


Slot Combine 


624 


2496 


936 


3432 



15 
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Diversity Selection Block 

The diversity selection block 138 selects one of the outputs of Antenna A and 
Antenna D (the two antennae locate at the edges) as one of the two inputs into the 
existing TRX. 



10 



15 



20 



AOA Estimation Block 

This block estimates the angle of arrival (AOA) of the strongest path that is 
used for the downlink beanifonning. The covariance matrix XX calculated in the 
beam former block is re-used in this block (this connection not shown in Fig. 10). But 
as there is only a 26 known bits sequence, preferably 4 samples per symbol data is 
used for this block Where X is given by: 

~£1(62) £2(62) £3(62) £4(62)~ 
£1(63) £2(63) £3(63) £4(63) 
£1(64) £2(64) £3(64) £4(64) 

_£1(87) £2(87) £3(87) £4(87)_ 
AndXX = X*.X 

Form a Hermitian Toeplitz matrix by using XX with the following procedures 
Zo = [XX(1 ,1) + XX(2,2) + XX(3,3) + XX(4,4)] / 4; 
Zj = [XX(1,2) + XX(2,3) + XX(3,4)] / 3; 
Zz = [XX(1 ,3) + XX(2,4)] 12; 
Z 3 =XX(1,4); 



ZZ = 



co«/(Z,) 



Z x 


z* 


z 3 ) 


z 0 


z, 


z 2 


conj(Zi) 


z 0 


z x 


conj(Z 2) 


conj(Z x ) 





(5) 
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• Do singular value decomposition of ZZ we may have ZZ = V A conj(V) T . 
where V is an orthogonal unit matrix formed by eigenvectors of ZZ and A is, a 
diagonal matrix formed by four eigenvalues. 

5 • Select the largest eigenvalue among the four and form the noise vector matrix 

by those eigenvectors not corresponding to the largest eigenvector. 

• Form a polynomial and find the root by looking up table or by decomposing 
the companion matrix. 

10 

• Convert the root into AOA in degrees and report it to Transmitter. 



Table 4 Complexity of AOA Estimation Block 





Real MPY 


Real Addition 


Total 


ZZ 


9 


6 




SVD 


256 


64 




Polynomial 
Roots 


2048 


10 24 




AOA 


1 


0 


3408 



Channel Estimation Block 

Multipath channels can be estimated by LMS method using the known 26 
training sequence. The four array outputs form four multipath channels, which contain 
all the information such as AOA, TOA, amplitude etc. Embodiments of the present 
invention fully exploit these multipath channels to achieve the best gain possible. We 
suppose the channel impulse has at least seven taps. Hence, we define a Toeplitz 
matrix S as 



15 



20 
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s(K 2 -6) s(K 2 -5) 



+ 6) 
s(K 2 ) 



to estimate seven time internals l/280-33K~ 3.7|is 



where s(Ki), s(Ki+l),..., s(K 2 ) are part of the known training sequence. Then each 
multipath channel impulse response can be obtained by solving the following linear 
equations: 



10 



15 



~cAl(6)~ 










, cAl(S) 








yKKi + 7) 


cAl(0)_ 




_N{K 2 -6)_ 






'ch2(6)~ 








~y2{K x +6) 


cA2(5) 


+ 


N{K X +1) 






cA2(0)_ 




N(K 2 -6) 




. > 2 (^ 2 ) . 


"cA3(6)~ 








>3(r i+ 6)" 


cA3(5) 


+ 






7) 


cA3(0)_ 




N(K 2 -6)_ 




. ^C^) _ 


"cA4(6)" 








>4(* 1 + 6)- 


cA4(5) 


+ 


N(K l+ i). 




^4(^+7) 


cfc4(0) 




A^C^ -6)_ 




. y^K 2 ) . 
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The explicit least mean square error solutions are: 



cM(6)' 
cM(5) 

cM(0) 

ch2(6) 
c/i2(5) 

cA2(0) 

cA3(6)" 
c/i3(5) 

ch3(0)_ 

cA4(6)' 
c/z4(5) 

c/z4(0) 



= (conj(SySy l conj(S) 



= (conj(S) T Sy l conjXSy 



= (conj(SySy l conj(S) 



6)' 

vlC^+7) 

. yl(^ 2 ) . 

>2(iT 1 +6)" 
^2(^+7) 

. J>2(* 2 ) 

'^3(^+6)" 
^3(r, + 7) 

>4( J RT 1 + 6) 
y4(^,+7) 



As the matrix S is formed by the known training sequence, the inverse matrix can be 
10 pre-calculated and stored in the memory. The total complexity for estimating the four 

multipath channels is given in Table 5. 



Table 5 Complexity of Channel Estimation Block 





Real MPY 


Real addition 


Total 


Chi 


4*(7*21 +7*7) 


2*(7*20+7*6) 




Ch2 


784 


364 




Ch3 


784 


364 




Ch4 


784 


364 


4592 



WO 2004/047333 



PCT/CA2003/001747 



26 

Referring to Fig. 13, there is illustrated in a block diagratn a transmitter for the 
Hybrid Space-Time Diversity System in accordance with an embodiment of the 
present invention. An existing transmitter 150 includes a BCCH TRX 152 and a 
plurality of channel transmitters 154. An RF/IF converter includes transmitter chains 
5 156 in the current deployed base station. The RF signal from 152 goes through the 

transmitter chain 156, which down converts the signal into base band signal, a 
watchdog function block 158 detects all the network information such as frame 
timing, training sequence, hoping sequence, which are fed into IF/RF block 118, 
antennas 160, and the receiver portion 104 of an applique intelligent antenna system 

10 (Fig 12). The plurality of transmitters 154 for communication channels are each 

having coupled to a transmitter chain 156 which down converts the signal and 
converts it to a digital signal for processing. A downlink beamfonner function block 
166 with weights from the receiver 168 processes the digital signals and provides an 
output to the IF/RF block 118 along with the BCCH output from the watchdog 

1 5 function block 158; 

A deployment of the Hybrid Space-Time Diversity antenna system is shown in 
Figs. 12 and 13, includes the following main components: 

20 L Three four-element antenna systems 160 (or, if deployed on a building in an 

urban setting, four four-element systems). 

• 120° scan angle for semi-urban or rural deployments 

• 90° scan angle for urban deployment 

• 60° scan or look angle for certain urban deployments 

25 

2. Certain electronic components on the tower or building, consisting of: 

• LNAs 

• Lightning arrestors 

• Converters to reduce the number of cables down the tower 
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3. Cabling to bring the RF down the tower or building; and DC power up the 
tower to feed the electronics 

4. User defined Shelter and base station. 

Further detail of an implementation of the antenna system include the 
following components: 



1 . Four antenna elements 
10 • Beam widths: 

o Choice of 120°, 90° or 60° in the azimuthally plane 
o 5.6° in the elevation plane 

• Physical size 

15 o Height: 6ft 

o Width: 14 inches 
o Depth: 6 inches 

• Gain: options; 

20 o Azimuths: 120°, 90° and 60°: 

o 17. 1 dB, 1 8.3 dB, and 20. 1 dB, respectively. 

2. Four LNAs and 1 st Mixers 

3. Four phase coherent receivers, producing I and Q outputs (Here they represent 
25 in phase and quadrature phase) 

4. Multiplexer 

TenXc Watch-Dog 158 

As the intelligent antenna system in accordance with an embodiment of the 
3 0 present invention may be hooked up to transmitters of various vendor' s TRX, the base 
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station information such as frame number, timing, timing advance, frequency hopping 
pattern may not be directly available. In this case a Watchdog function assists to get 
all this information when necessary. The Watch Dog function is assigned the 
following responsibilities. 

5 

1 . Regularly check slot 0 of BCCH carrier, particularly do a fast search for frame 
boundary by using both FB and SB. 

2. Decode SB to calculate Tl, T2, T3 ' and then derive exact frame number 

3 . Calculate the frequency-hopping pattern. 

10 4 - Decode BCCH norm or extension information, which might be carried in 

either, slot 0, or slot 2, or slot 4 or slot 6 depending on the deployed control 
channel combination. This information can be used to get Timing advance and 
therefore for downlink beamforming power control. 
5. Decode PCH. 

15 6 - Decode AGCH. This information along with the information acquired from 

RACH (initial AOA here) can be used for mobile positioning and therefore 
downlink beamforming. 
7. Decode NCH. 

20 Frequency Correction Channel 

Frequency Correction Channel (FCCH) is a downlink-broadcasting channel. It 
is carried by frequency CO (BCCH carrier) and always locates at slot 0. This burst is a 
constant burst with 0's fed into the whole slot. Therefore this burst causes a constant 
phase signal, in fact, the resulting signal is an unmodulated signal with a constant 
25 frequency C0(MHz) + 1625/24 (kHz). A mobile phone first refers to this frequency 

and adjusts its local oscillator (LO) to achieve a frequency synchronization with the 
BS. 
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This burst appears every 10 frames counting started with 51 frames cycle 
numerology. The Watch Dog performs a fast sliding correlation to obtain frame 
boundary information. 

Further detail on the FCCH Channel structure can be found in the GSM standard (ref. 
GSM 05.02) 



Synchronization Channel 

The synchronization channel (SCH) carries frame synchronization information 
and base station (BS) identification. After decoding this channel, a mobile terminal 
knows which BS connection to hook up and the exact frame number the BS is 
transmitting. The synchronization burst (SB) is always paired with the frame burst 
(FB) that appears just 8 slots later. In other words, it always appears at slot 0 of a 
frame next to the frame a FB appears. As the present intelligent system needs to 
decode this channel, we will detail this channel information format and channel 
structure in the following paragraphs! 

SCH Message Format and Bits Ordering 

The information carried in SCH is (a) the base station identity code (BSIC) of 
the base station, (b) Tl, T2, T3', three parts of the reduced TDMA frame number 
(RFN) as specified in TS GSM 05.02. The Figs. 14 and 15 illustrate the detailed 
arrangement this information in the message. Refer ETSI TS 04.08 for more detail. 

SCH Encoding 

The burst carrying the synchronization information on the downlink BCCH, 
the downlink CPBCCH for Compact, and in CTS the information of the CTSBCH-SB 
and the access request message of the CTSARCH, has a different structure. 
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it contains 25 information bits {d(0),d(l),..,, d(24)}, 10 parity bits {p(0),p(l),.„, p(9)} 
and 4 tail bits. 

The ten parity bits {p(0),p(l)„...,p(9)} are defined in such a way that in GF(2) 
5 the binary polynomial: 

d(0)D34 +...+ d(24)D10 + p(0)D9 +...+ p(9), when divided by: 
D10 + D8 + D6 .+ D5 + D4 + D2 + 1, yields a remainder equal to: 
D9+D8 + D7+D6+D5 + D4 + D3 +D2 + D+ 1. 

10 

Thus the encoded bits {u(0),u(l),...,u(38)} are: 
u(k) = d(k) for k = 0,1,...,24 
u(k) = p(k-25) for k = 25,26,...,34 
u(k) = 0 for k = 35,36,37,38 (tail bits) 

15 

The bits {e(0),e(l),..., e(77)} are obtained by the same convolution code of rate 1/2 as 

for TCH/FS, defined by the 

polynomials: 

20 G0 = l+D3+D4 

Gl = 1 + D + D3 + D4 

with 

e(2k) = u(k) + u(k-3) + u(k-4) ' 
25 e(2k+l) = u(k) + u(k-l)+u(k-3) + u(k-4) fork =0,1 77 ;u(k) = 0 fork<0 
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Synchronization Burst Transmission 





Bit Number 
(BN) 


Length 


Contents 

of field of field 


Definition 


5 


0-2 


3 


tail bits 


05.02 or below 




3-41 


39 


encrypted bits (e0..e38) 


05.03 




42 - 105 


64 


extended training sequence bit 


05.02 




106 - 144 


39 


encrypted bits (e39...e77) 


05.03 




145 - 147 


3 


tail bits 


000 


10 


148-156 


8.25 


guard period (bits) 


05.02 or below 



■ where the "tail bits" are defined as modulating bits with states as follows: 
(BN0, BN1, BN2) = (0, 0, 0) and (BN145, BN146, BN147) = (0, 0, 0) 

■ where the "extended training sequence bits" are defined as modulating bits 
15 with states as follows: 

(BN42, BN43 ... BN105) = (1, 0, 1, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 

0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 0, 

1, o, i, o, o, o; i, o, i, o, i, i, i, o, i, i, o, o, o, 

0,1,1,0,1,1) 

20 Frame Number Calculation 

After having decoded the SCH and mapped the bits into corresponding 
integers Tl, T2, T3\ then the frame number.FN can be calculated by 

FN = 51*((T3-T2) MOD 26) + T3 + 51*26*T1 
25 Where T3 = 1 0*T3 ' . For further detail see the GSM standard, Ref GSM 05 . 1 0. 
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Frequency Hopping Sequence Generation 

For a given set of parameters, the index to absolute radio frequency channel 
. number (ARFCN) within the mobile allocation (MAI from 0 to N-l, where MAI=0 
represents the lowest absolute radio frequency channel number (ARFCN) in the 
5 mobile allocation ARFCN is in the range 0 to 7 023 and the frequency value can be 

determined according to GSM 05.05 sec 2 with n= ARFCN), is obtained with the 
following algorithm: 

if HSN = 0 (cyclic hopping) then: 

MAI, integer (0 ... N-l) : MAI = (FN + MAIO) modulo 1ST 

10 else 

M, integer (0 ... 152) : M = T2 4- RNTABLE((HSN xor T1R) + T3) 

S, integer (0 ... N-l) : M' = M modulo (2 A NBIN) 

T = T3 modulo (2 A NBIN) 
ifM'<Nthen: 
15 S = M' 

else 

S = (M'+T) modulo N 

MAI, integer (0 ... N-l) : MAI = (S + MAIO) modulo N 

End 

20 NOTE: Due to the procedure used by the mobile for measurement reporting when 

DTX is used, the use of cyclic hopping where (N)mod 13 = 0 should be avoided, 
where: 

Tl R: time parameter Tl , reduced modulo 64 (6 bits) 
T3: time parameter, from 0 to 50 (6 bits) 
25 T2: time parameter, from 0 to 25 (5 bits) 

NBIN: number of bits required to represent N = INTEGER(log2(N)+l) 
A : raised to the power of 

xor: bit-wise exclusive or of 8 bit binary operands. 
MIAO : Mobile Allocation Offset Index (0 to N-l, 6 bits). 
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Downlink Beamformer for POC 

In order to simplify the implementation complexity, downlink beamfonner for 
POC will, be a fixed beam rather than adaptive one. Each sector has seven pre- 
designed fixed beams, a respective one pointing to -45, -30, -15, 0, 15, 30, 45 
5 degrees. The corresponding weight vectors are named as 

W a = [w a (l)w a (2)w a (3)w a (4)],W b = [w b (l)w b (2)w b (3)w b (4)],Wc- [w c (l) w c (2) 
w c (3)w c (4)],W d = [w d (l) w d (2) w d (3) w d (4)], W e - [w e (l) w e (2) w e (3) w e (4)], w; 
= Ml) Wf(2) Wf (3) Wf (4)], W g = [w g (l) w g (2) w g (3) w g (4)]. 

0 

Referring to Fig. 16, there is graphically illustrated downlink beam patterns. 



